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A note on the scale evolution of the ETQS function Tp(x,x) 
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Abstract 

We reexamine the scale dependence of the ETQS (Efremov-Teryaev-Qiu-Sterman) twist-3 matrix 



with the method 



C*"} ■ element which has been studied already by the four different groups with conflicting results lHJ]. 

We find that we can in fact reproduce the results of 



when we treat some 



r- 1 

subtleties with greater care, thus easing the mentioned conflict. 
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The ETQS matrix element plays an important role for the theoretical description of 
transverse single spin asymmetries (SSA) in the framework of collinear factorization. The 
control of Q-evolution is not only necessary to describe QCD dynamics correctly and to 
reduce the dependence of theory predictions on the factorization scale adopted, but such 
evolution equations give also insight into the functional form of higher-twist distribution 
functions. The idea there is to start evolution at a low scale and use the fact that the 
resulting form at a high scale is only little dependent on the low-scale input [5J. The latter 
is especially important in view of the limited experimental input of has to determine these 
functions. 

The corresponding calculation was done in Refs. \MiM- However, the result obtained in 



Ref. 



4] differ from that derived in Refs. l|-[3| by two extra terms. It was settled rather 



u 



easily that one of these terms is due to a Feynman diagram which was missed in Refs. |l|-|3]. 
The second additional term in [4| which is proportional to 5(1 — z) could not be reproduced 
by the other calculations so far. We show in this short contribution, how this term arises 
within the formalism of Ref. [2| due to a rather subtle fact related to the non-commutativity 
of a certain limit and a certain integration. We now hope to be able to do this calculation 
consistently in the light cone gauge. 

The ETQS function Tp is defined through the following matrix element, 

^r^e^ p+ y-(PS\Mo)9F^(yT)My-)\PS) = fr F (x,x)eTs ± j (i) 

In Ref |2j, the light-cone gauge (A + = 0) with the retarded boundary condition, i.e., 
Aj_(— oo~) = was chosen such that Tp(x,x) can be rewritten as, 

T ^) = ( P^je^ p+ y- (PS\m^TS±uid^ a (y-)\PS) . (2) 



Itt 2 M 

To calculate the splitting function, one has to take into account the contributions from 
the operators (^d±_ip) and (ipA±ip^, because they are of the same twist. We plot the 
Feynman diagrams contributions for the real gluon radiation in Fig. 1, where (a) is the 
contribution from the partial derivative acting on the quark field, and (b — d) are those from 
A± contributions. Virtual corrections only contribute to the contribution proportional to 
NJd±ip) . Their contribution is the same as for the quark self energy correction. 

Following the procedure presented in Ref. |2[, we perform a collinear expansion for the 
hard scattering part to calculate the contribution from Fig. 1(a). The linear k±_ expan- 
sion term combining with the quark field will lead to the quark-gluon correlation function 
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FIG. 1: Real gluon radiation contribution to the evolution equation for the ETQS function Tp{x, x). 
Crosses in fig. (a) and horizontal bar in fig.(b) indicate k± flow and special propagator, respectively. 

Tp(x, x). In the collinear expansion in terms of k±, we can fix the transverse momentum of 
the probing quark (l q ) or the radiated gluon (l g ), because of momentum conservation and we 
are integrating over them to obtain Tp(x, x). We have also checked that they will generate 
the same result. Following the Ref. 2], we fix l g in the collinear expansion to simplify the 
calculation. 

For the A± contribution, we notice that F +fl = d + A^_ in the light cone gauge. Therefore, 
one can relate the corresponding soft matrix to the correlation function Tp{x,X\) in the 
following way, 

C rim rlni 

(PS\Mo-)^s ±u gF+^)My-)\PS) 



i i dy_dyj_ e i X1 p+ y - pi ( x - Xl )p+y- 

x — x\ + ie .1 Att 



dy dy 1 

4ir 



P+e^ p+ y-e^- x ^ p+ ^(PS\Mu~)^7s±vgA ± ^)My-)\PS) . 



(3) 



In above formula, the soft gluon pole appearing in the first line is generated by partial 
integration. The pole prescription has been determined by our choice of a retarded boundary 
condition. For the same reason, we have to regulate the light cone propagator in a consistent 
manner. The gluon propagator appearing in Fig. 1(c) in the light cone gauge with the 
retarded boundary condition is given by, 



D a P(l) 



a 6 



l a nP + n a P 



(4) 



I 2 + it \ I ■ n + ie 

where I is the gluon propagator momentum flowing out from the quark-gluon vertex in 

Fig.l(c). 



We now deviate from the original calculation [2] in two ways: 



i) in [2] the integral f , dx. 



x'5(x') 



x- 



was simply neglected; 



ii) one of the two absorptive parts of the free propagator was not taken into account. 

We will discuss next these two points in more details, arguing that the neglected contri- 
butions have to be taken into account. When computing the hard pole contribution from 
Fig. 1(c), for the left cut diagram, one has 

/L> r ' r (x + x B )(2x g -x g )~ 



o-i(l)|hp-left / \ 

J- IT Il7;_ 1 f„"\ [•LB I 



«s f dx f dl gL C A 
4tt.L x .L.,* 9 ! 2 2 



2x 2 g 



F iFig.l(c) 

-" jx b ~ ux' g ~g± - L —" g J 

where x' g = l-n/p + with x' g = xb — x + x g . By noticing that J x , dx g Xg Jt 9 
rather than zero, and summing left and right cut diagrams, one obtains, 

Os_ r dx dl 2 g± C A 
2vr J XB x P g± 2 



5(x')T F (x - x g ,x)(5) 



S(x B 



T (i),h P 

1 F I Fig. 1(c) 



{xb) 



l + Z 



T F (x B , x) - 5(1 - z)T F (x B , x b ) 



x 



(6) 



where z = xb/x. The second term is missing in Ref. ^ 

Next we discuss the second contribution which was overlooked in Ref. |2|. Since we work 
in the light-cone gauge with retarded boundary condition, the free propagator possesses two 
absorptive parts [6j, 



discD a \l g 



2k9(IW q ) 



L g±- 



^e(i° q )5(it) { -^— 2 — ii (?) 
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In [2j only the first absorptive part was taken into account. In order to carry out the 
calculation in a consistent manner, one must include the contribution from the second part. 
However, if one still picks up the same imaginary part as we did above, this contribution 
will cancel between the different cut diagrams as it happens when both gluon lines go on 
shell. On the other side, the additional imaginary part may come from the artificial pole 
which appears in Eq.(3). Such pole-absorptive part combination gives the contribution, 



T. 
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L g± J*b 
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dl; 



dx g 5(x g — x')S(x g 



x 



C A 



2(2x B - x' g )(x g + x' g ) 
2{2l-x B + l 2 g± )x' g " (21-xb + ZJJ 



2x B l 2 g± 



T F (x — x g ,x) (8) 



Integrating over x g , I and summing the two cut diagrams, we obtain, 



T 
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dx Ca 
~^2 



dy- 



S(l-z)T F {x B ,x B ) (9) 



Taking into account the contribution from the second part of the absorptive part, Eq.(20) 
in the Ref. [2J should be modified as follows, 
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T F (x,x) 



(10) 



Collecting all pieces, we eventually arrive at the following scale evolution equation for 
T F (x,x), 



dT F (x B ,XB,pt> 2 



91n/i 2 



2tt 



■I'D 



dx 
x 



C h 



+ ■ 



C A [l + z 



T F (xz, x) 



\-z 
l + z 2 



-5(1 - z) }T F (x,x) 



T F (x, x) — 25(1 — z)T F (x } x) + T F (xz } x) 



(11) 



which coincides with the result given in Ref. |4]. 

As shown above, the missing boundary term —25(1 — z)T F (x,x) appears to be a generic 
problem which might have far reaching consequences. In principle, all previous calculations 
involving hard gluon pole contributions might need to be reexamined. However, one can ex- 
pect that the observed matching between the TMD factorization and collinear factorization 
at intermediate transverse momentum will not be affected by this extra term. 
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